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Abstract 

In this work we examine the effect of phase-space noncommutativity on some typically quantum 
properties such as quantum beating, quantum information, and decoherence. To exemplify these 
issues we consider the two-dimensional noncommutative quantum harmonic oscillator whose com- 
ponents behavior we monitor in time. This procedure allows to determine how the noncommutative 
parameters are related to the missing information quantified by the linear quantum entropy and 
by the mutual information between the relevant Hilbert space coordinates. Particular issues con- 
cerning the thermodynamic limit of some relevant properties are also discussed in order to evince 
the effects of noncommutativity. Finally, through an analogy with the Zeeman effect, we identify 
how some aspects of the axial symmetry of the problem suggest the possibility of decoupling the 
noncommutative quantum perturbations from unperturbed commutative well-known solutions. 
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I. INTRODUCTION 



Loss of quantum coherence, quantum entanglement, quantum beating, wave function col- 
lapse and several related issues are at the heart of quantum features that must be necessarily 
considered in order to construct a suitable quantum measurement theory [IH3]. 

These effects allow for a better understanding of the frontiers between quantum and 
classical descriptions j2] as, for instance, about the generation of superposition states of 
mesoscopic systems, known as Schrodinger cat states [HE]. 

The central question examined in this work concerns typically quantum mechanical fea- 
tures in the context of a phase-space noncommutative (NC) extension of quantum mechanics 
(QM) [6H9]. More particularly, the effect of phase-space noncommutativity on phenomena 
such as quantum beating, missing information and decoherence. 

To start with, noncommutativity in the configuration space has been firstly suggested 
as a way to regularize quantum field theory [10]. Furthermore, it is a long-held belief 
that in quantum theories which include gravity, the nature of the space-time is modified 
at the Planck scale, and that noncommutativity must be considered at high energy scales. 
Noncommutative geometry is also encountered in the context of string theory /M-theory [TIT 
[13] . We also mention that several aspects of NC theories have been explored along decades, 
including noncommutative gauge and field theories [T^ - fl9] . 

In the quantum mechanical context, NC extensions have been discussed with respect to 
the quantum Hall effect [2"0l |2"T] . for electrons in a magnetic field projected to the lowest 
level in the Landau problem in the phase-space [2], for the 2-dimensional (2D) quantum 
harmonic oscillator in Cartesian coordinates [22], and in the gravitational quantum well 
for ultra-cold neutrons [231 121] . The NC QM is supported by canonical extensions of the 
Heisenberg-Weyl algebra [231428] . The theory lives in a 2<i-dimensional phase-space where 
time is required to be a commutative parameter, and coordinate and momentum variables 



obey a NC algebra (c. f. Eq. (gT))). 

It is important to stress that phase-space NC extensions of QM have striking features. 
These include putative violations of the Robertson-Schrodinger uncertainty relation [8], a 
feature that has somehow been observed in a recent experiment [29] . Moreover, the regular- 
izing features of phase-space NC QM are remarkable and shown in minisuperspace quantum 
cosmology models [30J, and in black-hole physics where a L 2 wave-function can be found 
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and the singularity problem can be resolved [3TJ [32] . 

The phase-space NC extension of QM is more suitably formulated in the Weyl-Wigner- 
Groenewold-Moyal (WWGM) formalism for QM [3"3Tf3"o*] . Its framework based on the 
WWGM understanding of the QM statistical structure has been developed in the last years 
[6j [36VH7] . Thus, it is natural to use the WWGM formalism in our study of the 2D quan- 
tum harmonic oscillator. As will be seen, NC introduces several conceptual subtleties in 
describing the dynamics of quantum systems and affect some typical quantum phenomena 
such as quantum beating, ensued missing information and wave function distortions. The 
quantum beating related to the interchange of information between the Hilbert sub-spaces 
of the problem are, as will be seen, exclusively due to the NC properties. Not less important, 
the missing information due to probability distribution distortions may be evinced through 
the coupling of NC degrees of freedom to a thermal reservoir. 

The outline of this manuscript is as follows. In section II, we review the Weyl-Wigner 
formalism and its applicability in describing the NC extension of the QM in the phase-space. 
We report about the formal procedure for obtaining the energy spectrum and the Wigner 
storgenfunctions of a NC quantum system through the Groenewold-Moyal star-product in 
the NC phase-space that exhibit both coordinate and momentum noncommutativity. In 
section III, we apply the formalism to calculate the NC properties of the quantum harmonic 
oscillator on the plane. Besides reproducing some preliminary results on the stargenvalue 
problem of the harmonic oscillator in the NC phase-space, our emphasis is on obtaining 
the phase-space time-evolution that leads to the time-dependence of the Wigner function. 
This provides a fruitful framework for investigating NC effects. We show how noncommuta- 
tivity potentially affects the quantum coherence of a state vector that is obtained through 
an arbitrary quantum preparation procedure. Hints about the possibility of observing non- 
commutativity through quantum beating effects are discussed. The missing information 
quantified by the linear quantum entropy is discussed in section IV. The quantum entropy 
destroyed by the NC character is quantified and a well-defined prescription for the wave 
function distortion in thermodynamic equilibrium is obtained. Results for the thermody- 
namic limit of the internal energy, Boltzmannn entropy and heat capacity are derived. The 
quantum issue of the 2D NC harmonic oscillator is reformulated in section V. It is shown 
that some aspects of the axial symmetry of the problem allow for a factorization of the NC 
effects. Obtaining a well-defined quantum number and a factorable wave-function compo- 
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nent for noncommutativity are shown to be relevant for quantifying the NC effects on the 
2D quantum harmonic oscillator. Finally, we draw our conclusions in section VI and address 
some related issues concerning the anisotropic NC harmonic oscillator, and possible extra 
linear contributions to the potential in space coordinates. 



II. THE WEYL TRANSFORM AND THE WIGNER FUNCTION IN THE 
PHASE-SPACE 

The predictive power of QM relies on its formalism to connect quantum observables to 
expectation values. In this respect, the Wigner function and the Weyl transform establish 
an alternative formulation to better known Schrodinger and Heisenberg frameworks. This 
formalism provides an interesting insight into the nature of QM and its classical limit. The 
Wigner representation has also the operational advantage of exhibiting all the information 
about the state vector in the phase-space. This contrasts with the conventional QM approach 
which provides information about position, or momentum, but not both together. This 
formulation is akin to the formalism of statistical mechanics. 

The inception of the Wigner function is the definition of the Weyl transform of an operator 
0(q,p) through the expression 

O w (q,p) =Jdy exp [ipy/H] (q-y/2\6\q+y/2) = Jdy exp [-iqy/h] (p-y/2\0\p+y/2), (1) 

in both coordinate and momentum basis, where operators, q and p, are converted into c- 
numbers, q and p. The fundamental property of the Weyl transform is that the trace of the 
product of two operators, 0\ and O2, is given by the phase-space integral of the product of 
their Weyl tranforms [351 EH] as 



Tr 



{<?.p} 



0^2 



h- 1 dqdpOT(q,p)0 2 V (q,p), (2) 



where h is Planck's constant. The trace properties can be related to the density matrix 
properties of the state vector described by p: 

Tr iq>p} [ P 0] = (O) = h- 1 J Jdqdpp w (q jP )O w (q,p). (3) 

By identifying the density matrix p with |^)(^|, one can define the Wigner function, 

W(q,p) = h- l p w = jdy exp [ipy/h] V(q - y/2) V*(q + y/2), (4) 



which can be naturally generalized to a statistical mixing, such that the expectation value 
of an observable O can be computed through 

(0} = J JdqdpW( q ,p)0 W (q,p). (5) 

By identifying *&(q) with 

^(q) = h~ x Jdp exp [zpg/ft] (6) 

one straightforwardly obtains probability distributions for q and p as 

fdqW(q,p) = $*(p)$(p) and JdpW(q,p) = (7) 

i. e. the Wigner function represents the distribution in the phase-space represented by ^(q), 
and can also be computed from as 

W(g,p) = Jdy exp[-iqy/H)${p-y/2)<$>*{p + y/2). (8) 

Additional properties related to the density matrix theory can be obtained from the above 
prescription [49, 50J. Our interest is mainly concerned with the definition of Tr[p 2 } as 

Tr[ P 2 ] =h~ 1 J jdqdpW(q,p) 2 , (9) 

that demands for an extra normalization factor, h~ l , in order to keep consistence with the 
density matrix theory that sets Tr[p 2 ] = Tr[p] = 1 for pure states. 

Providing a distribution in the phase-space and a representation for the expectation val- 
ues of quantum observables is not, nevertheless, an exclusive property of Wigner functions. 
A similar procedure could be implemented through other function candidates to provide 
enough information about the quantum systems. Given that it admits positive and nega- 
tive amplitude values, the Wigner function cannot be strictly interpreted as a probability 
distribution. On its hand, the Husimi distribution, for instance, is an alternative which 
ensures non- negative values and admits a probability interpretation [ST]. This distribution 
is largely used in quantum optics [53] and to study quantum effects in superconductors [54J . 
Likewise, to describe light in optical phase-space, a third proposal, the Sudarshan-Glauber 
P representation [55] is a common tool. 

Despite these examples, due to its simpler formulation, the Weyl- Wigner formalism is the 
most often procedure used to overcome some obstacles in quantum dynamics. To provide 



an example of its functionality, let us recover the commutative Heisenberg-Weyl algebra of 
the ordinary QM, 



o. 



Qii Pj 



ih5; 



I.J- 



p p 

1 n 1 3 



0. 



i,j = l,...,d . 



(10) 



The Weyl transform can be extended to a d-dimensional phase-space, {Q, P}, such that any 
quantum operator, 0(Q, P; t), in terms of c-numbers, can be written as 



O w (Q,P;t) 



dx. dyF(x, y ; t) exp 



h 



(x-P + y-Q) 



where, through the d-dimensional generalization of the Weyl transform, one has 



F(x,y;t) = h- d Tr {Q , P} 



0(Q,P;t) exp 



h 



(x-P + y-Q) 



(11) 



(12) 



The algebra in Eq. (10) is recovered through the introduction of the Groenewold-Moyal 



star-product defined on the space of commutative functions as 

"A 



h-k g = exp 



h(r)g(s)\, 



(13) 



where A is a suitable constant. 

The main feature of WWGM phase-space formalism is that 



;*|0 1 (Q,P;*)0 2 (Q,P;*)|*) = J jdP dQp w (Q,P;t)0^(Q,P;t) *0? (Q,P;t), (14) 



where 



p w (Q,P;t) = h- d I dzexp[^z-P](Q-^|p|Q + ^), 



(15) 



which shares the same properties of the Wigner quasi-probability distribution function from 



Eq. (11), once the density matrix for a pure quantum state, p = \^}(^/\ is introduced. The 
Moyal bidifferential star-operator is given by 

iH , 



exp 



such that the dynamics of the Heisenberg operator is denoted by 

O w (Q, P; t) = -- \O w (0), H W (Q, P; f)l = Uh w ★ O w - O w * H w ) 
n n 



(16) 



(17) 



which follows naturally from the WWGM formalism. From the above construction, the 
commutative variables P and Q satisfy the Hamilton equations of motion, namely 



P 

Q 



h 

i 
h 



[p w (0),H w (Q,P;t)] ic = -V Q H W ; 
[Q w (0),H w (Q,P;t)]^+V P H w , 



(18) 



and, therefore, the corresponding classical dynamics interpretation is inherent from the 
formalism. 

Once again, the Wigner function projection on configuration and momentum spaces yields 
the quantum mechanical configuration and momentum probabilities, respectively. Although 
somewhat beyond the scope of this work, several additional properties can be recovered from 
the WWGM formalism |6j [221 HE] • For instance, by making use of the integral representation, 
one has 

Or(Q,P)*Of(Q,P) = (2nh)- 2d [ ... fdP' dP"dQ' dQ"OY(P',Q')0^(P",Q") 

J J (19) 
exp[--(P ■ (Q' - Q") + P' • (Q" - Q) + P" • (Q - Q'))], 

and obtains that 

dP dQ H w (Q,P)*p w (Q,P) = [ [ dP dQ H w (Q,P)p w (Q,P) = E, (20) 




which shall be useful in the following siargenvalue problems involving NC variables. 
Turning to the extended NC algebra, the following commutation relations are satisfied 

[qiAj] = iOij, [VhPj] = ittiij, \PhPj] = iVij, i,j = 1, -,d (21) 

where r/y and 9{j are invertible antisymmetric real constant (d x d) matrices, one can define 
an equally invertible matrix 

Sjj = 5ij + -^9ikr]kj, (22) 
if for any matrix elements, rj and 6, their product is smaller than h 2 [6], 

Or, < h 2 . (23) 

Under a linear transformation this algebra can be mapped into the usual Heisenberg-Weyl 



algebra as in Eq. (10) via the Seiberg-Witten (SW) map [13], which can be cast in the form 

ft = A, j Qj + Bijtlj ^ = dj Qj + Dijflj, (24) 
where A, B, C, D are real constant matrices. The above transformation constrained by the 



NC relations from Eq. (10) is easily shown to obey the following matrix equations 



AD T — BC T = Idxd AB T -BA T = ^e CD T - DC T = , (25) 

n n 
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where A, B, C, D, ©, N are matrices with entries Ay, B^, C^, Dij,8ij, rjij and the superscript 
T stands for matrix transposition, respectively. 

The linear transformation ensures that the abovementioned NC algebra admits a repre- 
sentation in terms of the Hilbert space of ordinary QM. In the noncommutativity context, 
the Wigner-function is much more than a visual insight into the meaning of quantum states. 
As it will be pointed out in the following section, the Weyl-Wigner formalism is actually the 
most natural approach for applying the NC formalism. 

III. THE TIME-DEPENDENCE OF THE 2D NC HARMONIC OSCILLATOR: 
QUANTUM BEATING AND DECOHERENCE 

The generalization the Weyl functions and the Wigner distributions to the NC phase- 
space and the general features of the Weyl-Wigner map, which includes its covariant general- 
ization to d-dimensional systems, have been formally developed in Ref. [6J. The generalized 
siar-product, the extended Moyal bracket and the NC Wigner function have been mathe- 
matically formulated, and it is found that the extended Weyl-Wigner map is independent 
of any particular choice of the SW map. 

It is not our purpose to re-discuss here the mathematical grounds of the formalism and 
its general applicability to quantum systems, as performed in [6], but in order to relate the 
noncommutativity properties with observable quantum phenomena, it will suffice to apply 
the results to the well-known harmonic oscillator on the plane. 

The quantum mechanical problem of the 2D harmonic oscillator is the most accessible 
system for which the NC phase-space properties can be probed [SI 122], and quantum effects 
quantified. 

Thus, let us consider the quantum Hamiltonian, 



on the NC x — y plane, with spatial and momentum noncommutativity being set by 



##o(q,p) 




2 2 

H — muj q 



(26) 



[ft, Qj] = iOe 



[QuPi] = ibjft, \Pi,Pj) = iriUj , hj = 1,2 



(27) 
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with €ij = —€ij, such that the map to the commutative operators is given by 

a ( QvY 1/2 ( & 
g < = H J [ q * + 2\fIh eijPj 

n^A(l-|) (A -^^), (28) 

through the SW map, 

9 7] 

q% = Aft - — eijUj , pi = fiUi + ^-^eijQj , (29) 
which is invertible when the parameters A and /i are constrained by the relationship 

^ = AMl-A/i), (30) 
with 9rj < h 2 , and with the corresponding Jacobian reading 

||| = (de t ^ = 1 -|. (31 ) 

The Hamiltonian in terms of the commutative variables, Qj and ITj, reads: 

2 

tf]f (Q, n) =a 2 Q 2 + P 2 n 2 < u \\,Q r (32) 



where 

a" = — h 



2 X 2 mu 2 i ?7 2 



^ a 2 mco 2 9 2 

8 = — — l 

' 2m 8A 2 ft 2 ' 

1= 2h mU+ 2mh' (33) 



Giving that, from Eq. (18), the commutative variables, Q and II, satisfy the Hamilton 
equations of motion, one obtains the following set of coupled first-order differential equations, 

i 



life = -j. [n fc , H™ ] = -2a 2 Q k --ye jh ILi, 

Qk = -\[Qk, H% ] = 2f3 2 Ii k - 1 e ]k Q J , k,j = l,2 (34) 

so that Q and II may be interpreted as classical dynamical variables within the Wigner 
formalism. The above equations can be re-written as two uncoupled forth-order differential 
equations as 

n* = -2( 7 2 + 4a 2 /3 2 ) ii fc + ( 7 2 - 4a 2 /3 2 ) IT fc , 

Q k = -2(7 2 + 4a 2 /? 2 ) Q k + ( 7 2 - 4a 2 /3 2 ) Q k , (35) 



from which one gets the solutions 

Qi(t) = x cos(fit) cosHt) + y cos(£lt) sin(jt) + — [ir v sin($7t) sin^i) + ir x sin(Ot) cos(7i)l , 

a 

Q2(t) = y cos(Oi) cos(7t) — x cos(fii) sin(7f) — — [ir x sin(Oi) sin(7i) — n y sin(Ot) cos(7i)] , 

OL 

IIi (i) = t^x cos(rit) cos(7t) + ir y cos(f2i) sin(7i) — — [y sin(Qi) sin(7t) + x sm(Qt) cosfat)] , 

P 

OL 

Il2(t) = % cos(J7t) cos(7t) — 7r x cos(fii) sin(7i) + — [x sin(Oi) sin(7i) — y sin(Ot) cos(7t)] , (36) 
where x, y, ir x , and 7r y are arbitrary parameters, and 

Q = 2afi = uy/(2Xfi-l) +e 2 , (37) 

with 



1 ' " (38) 



2h 



mujQ H 



moo J 



so that 7 = use. By conveniently setting A = /i = 1, one sees that Q ~ w(l + (9(e 2 )), so 
that the NC parameters, 9 and rj, introduce second-order modifications onto Q. Likewise, 
the modifications due to 7 = ue correspond typically to first order effects. Notice that by 
setting e = one recovers the solutions for the 2D harmonic oscillator with uncoupled x — y 
coordinates. The above results lead to time-invariant quantities expressed by, 
2 



i=l 

2 

^2 ( e ijQi(t) n?(0) = xn y -yii x , (39) 

i=l 

which will be useful in the subsequent discussion. 

The time evolution of the phase-space coordinates, ( Q\ (t), Hi(t)) and (Q2{t), n 2 (i)), for 
the NC harmonic oscillator are depicted in Fig. [TJ For convenience, we define an auxiliary 
variable, e, as 

e = 2 = y £ » (40) 
in order to perform a non-perturbative analysis of the results. If e can be written as a 
rational number, one encounters a beating effect along the phase-space trajectories. The 
same does not occur for non-rational values of e (c. f. the last plot in Fig. [I]). 

Following the notation of Section II, the stargenfunctions for the Hamiltonian problem 



from Eq. (32) are obtained from the stargenvalue equation, 

a no * pZu 2 (Q, n) = E nun2 P ^ n2 (Q, n), (4i) 
10 



which, from the analysis developed in Ref. [22J, results into 
f—1 y«i+«2 



cxp 



i(> 2 + - n2 



L ni (n + /h)L n2 (n„/h), (42) 



where L° are the associated Laguerre polynomials, ri\ and n 2 are non negative integers, and 

2 

■ 1 ..,./, , . 



n± = ^ 2 + -n 2 T 2V( eiJ n t 



a 



(43) 



such that the energy spectrum is given by 



and one has 



Enuna = ^ [2a/3(™i + n 2 + 1) + 7(^1 - "2)] 



dQi / dnx / dft / rfn 2 P ^ jn2 (Q,n) = 1. 



(44) 



(45) 



From the results of Ref. [6], one realizes, using the inverse SW map Eq. (28 ), the constraint 
Eq. (30), and the normalization (detO) -1 / 2 = (1 — Orj/h 2 ) \ that the reported NC Wigner 



function is independent of A,//. The implications of the above results for the wave function 
time evolution are discussed in the following. 

A. Quantum beating and decoherence 



The dynamical evolution of a state vector described by a generic Wigner function, 
p w (Q, II), follows the motion of coordinates in the phase-space, {ft, IT}. One has only 
to ensure that each point of the Wigner function moves in the correlated paths, 1 -H- 2, 
as depicted for instance in Fig. [T] for the NC harmonic oscillator. This is reflected by a 
characteristic invariance property of Wigner functions that sets 



p w (Q,U)^p w ( Q, U, t) = p w ({x, 7t x }, {y, tt J; 0) 



(46) 



as one can obtain from the inverse transformation derived from Eq. (36) [33]. If the state 



vector preparation follows the prescription of the quantum numbers n\ and n 2 , from Eq. (39 ), 



one can verify that the condition prescribed by Eq. (46) results into stationary Wigner 
functions, p^ n2 (Q,n). 

However, the properties of the Wigner function are not constrained, in general, by the 
dynamics of the system. For the NC harmonic oscillator, an analysis using the properties of 
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the Laguerre polynomials (c. f. Eq. (42)) shows that some features are associated with the 
preparation of states rather than to the physics of its time evolution. 

One should notice that, for instance, by re-writing the quantum numbers for the 
storgenvalues through the constraint n = n\ + 712, the state vector can be prepared as a 
normalized statistical mixing with a well-defined quantum number, n, as 

(-11 " 



pf(Q,n) 



l + n 
1 f 



ni=0 
-l) n 



exp 



n \p a 



(n + /h) L° n _ ni (n_/ft) ,(47) 

ni=0 



7T 2 /l 2 l+n 

from which, by noticing that 

n 

E ^ (°) ^ ( & ) = L n d+l (« + & ) > 

711=0 

one obtains 



where 



a ~2 , ^2 



p a 



a 



(48) 



(49) 



(50) 



Results from section V suggest that such a normalized quantum superposition, p£ (Q, II), 
can be related to the solution of the commutative 2D harmonic oscillator written in terms 
of axially symmetric coordinates. By using the same Laguerre polynomial identity from 



Eq. (48), one can find the usual decomposition of p^(Q,Tl) into orthogonal Cartesian 
coordinates, x and y, such that 

(-l) n " 



PniQi 11 ) = Pn(() 



1 + Tl 



^ ] Pn x (£3 



P 



xl rji- tlx vsj/ 



(51) 



^=0 



represents the quantum superposition of normalized uncoupled pairs of one-dimensional (ID) 
harmonic oscillators, where 

1 exp[-e x , y /h]Ll tny (2£jh), (52) 



pZJC 



with n, 



n — n x , and 



S.r 



Oi 
P 



x H — 71 
a ■ 



a 



P 



(53) 
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where we have used Eqs. (46) and (50). 



Since the state vector described by Eq. (49) does not exhibit any NC effect from Hamil- 



tonian Eq. (32), then a quantum system can be arranged in order to exhibit the properties 
from the state described by Eq. (49). However, even for the NC phase-space, p^ (£ 2 ) de- 



scribes a stationary state vector (c. f. Eq. (39)). The same is not true for the composition 
of state vectors defined by, 



pZ, ny ( Q, n; t) = pZ (^) pi (g) -+ p W(Q } ri; t)p™(Q,Il;t). 



(54) 



For the NC scenario described through the dynamics arising from Eq. (36) constrained by 



the invariance conditions, Eqs. (46) and (50), one can obtain the explicit time-dependence 



for the state vector from Eq. (54) by noticing that the time evolution of i^xy) m ^ x ^ ne 
components from x and y phase-spaces since one has that 



t ~ 



C 2 



a 




+ 


i 


t 


11 


a 


a 




+ 


P 


t 




a 



t + ^ 

a 

t + -n| 

a 



a 3 

-QiQ 2 + -U 1 U 2 )sin(27t) 

p a 



a 



.(55) 



The above defined state vectors are stationary for the commutative case (7 = 0), and they 
have indeed a time dependence mediated by 7. 

To show how the noncommutativity affects the state vectors prepared as pure states 
described by p^ x n , one looks over the time evolution of the traced-out Wigner functions 
(density matrices) in the corresponding phase-space. The traced-out Wigner functions are 
defined as 



p^iQ^t) = Tr {2} [^(ftnjt) 

/ + °o +00 
dQ 2 / dli 2 pZ^ ny (Q,U;t), 
00 — 00 

pW ny (Q 2l Il 2 ;t) = Tr {1} [p^ ny (Q,U;t)_ 

to „ + oo 

dQ 1 / cflli P Z, n AQ,n;t). 



(56) 



(57) 



i. e. tracing out over Q 2) i and 112,1 means integrating p^ x ny over these variables, so that the 
resulting Wigner function in the Q\ >2 — Hi j2 plane is obtained. Assuming that we look over 
the time evolution of the state vectors along the x direction, we consider the corresponding 
Wi gner functions, pn},n y {Qi,^-i'-,t) as shown in Figs. [2] and |3j At time r = 0, the Wigner 
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function corresponds to a description of the ID harmonic oscillator stationary state with 
quantum number n x . The quantum beating effect due to the NC parameter 7 is obtained 
by looking over time intervals multiples of 7r(8eo;) _1 . The NC parameter e = 7/Q modifies 
the usual time behavior of the commutative harmonic oscillator by generating the quantum 
beat with frequency Ub ea t = 27 = 2eu. Notice that the initial quantum state in the Q\ — TEi 
phase subspace collapses into the traced-out quantum state, by reproducing its wave pattern 
at r = Tr/(2eu). The results depicted in Figs. [2] and [3] do not depend quantitatively on the 
auxiliary parameter, e, since we have chosen scale independent values for r. 

To quantify the effect of the NC parameter, 7, over the time evolution of the state vectors 



like those from Eq. (54), one computes the linear entropy defined by 
S 1 (t) = 1 - ^Tr {1} \(Tt {2} \pZ(Q,n;t) 



s 2 (t) 



1 - ~H Tr ^ 
2tt 

1 " T Tr « 

2tt 

1 - ~^Tr {2} 



S l2 (t) 



h 

2tt. 
47rJ 



-Tr 



{2} 



Tr {1} 



p(2) 

' Tlx i^ly 



o w ( 



Q,n,t) 



Tr {l} 



Tr 



{2} 



(<?2,n 2 ;t) 



pZ >ny (Q,U;t) 



through which the quantum mutual information is defined by 



Ii 2 (t) = S 1 (t) + S 2 {t) - S 12 (t) = I 21 (t), 



(58) 



(59) 



which quantifies the correlation between x(<H- 1) and y(<H- 2) states. The mutual information 
is a measure of the correlation between subsystems of a quantum state. If two variables 
x and y are assumed to be uncorrelated, the quantum mutual information measures the 
discrepancy in the uncertainty resulting from this possibly erroneous assumption. By setting 
7 = 0, the mutual information, Ii 2 (t), as well as all the above defined entropies vanish, once 
Pr^,ny ( Qi t) reproduces the product of two uncorrelated pure states. If the NC effects 
vanish, the state vector evolves like a stationary pure state. 

Figs. [2] and [3] also exhibit a swapping dynamics due to NC effects. It works as if one had 
selected a specific coupling between oscillators on the NC plane, which would establish the 
dynamics of interchange of states or state-swap between them. l\ 2 quantifies such a peculiar 
behavior as found, for instance, in a weakly coupled QED cavity [56] . 
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In addition, if a time scale, r, is not sufficiently large to close the phase-space orb- 
ifolds (c. f. Fig. Q) for Qi(r) and Ilj(r), the quantum beating can be reinterpreted as a 
scale-dependent decoherence effect as qualitatively illustrated by Fig. |5j Depending on the 
magnitude of the NC parameter, 7, the beating frequency gives rise to a decoherence time, 

Tcoh = 7T/7. 

The relevance of these results in terms of its experimental feasibility/detectability may 
depend essentially on the parameter 77. By setting rj = m 2 co 2 6 one has hy = rj/m and 
the coherence time reduces to r co ^ = mh/2i], such that only for small values of the mass, 
m, NC can play a relevant role in evincing decoherence / quantum beating effects. The 
corresponding coherence length related to r co h can be defined as 

pi hi hw hm h 2 u 

coh m f Am 27 2-Km 2rj 77 ' 
where it has been assumed that c = 1. 

IV. MISSING INFORMATION AND THE THERMODYNAMIC LIMIT 

A further issue to be considered in this work refers to the effect of noncommutativity on 
the loss of quantum coherence, the ensued missing quantum information and its imprint on 
the thermodynamic limit of quantum states when variables like internal energy, Boltzmannn 
entropy and heat capacity are investigated. 

When one looks over a state vector in the NC plane x — y, in comparison with solutions of 
the commutative problem, the missing information can be quantified through the quantum 
entropy of the entire system in the Hilbert space. The most useful variable for this analysis 
is the linear entropy, Eqs. (58), computed over a thermalized statistical mixing of m — n 2 
states: 

pY h (a;r u ki) = ^ 2 — — exp[-a] exp [-{a^ + a 2 n 2 )] pZ, n2 { r i, k) (61) 

711,712=0 

with 

a s = a{\ - (-l) s e), s = 1, 2, (62) 



where p^ in2 are the normalized state vectors given by Eq. (42), a = hVL/ (hsT), and N(a) 
is the normalization factor that contains the information from the corresponding statistical 
weights. 
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This variable corresponds to the state operator of the infinite multiple modes n\ and n 2 
of the 2D harmonic oscillator, in equilibrium with a thermal bath through an interaction 
Hamiltonian factorizable as Hi <8> %2 ® Hbath m the Hilbert space representation. 



One notices that the state Eq. (]61|) is quite suitable for analytical manipulations. By 

T 



introducing the coordinate parametrization, 




Qi and ki 




not 



(63) 



such that dr\ dki = dQi dUi, with i — 1, 2, and re-defining 



i=i 

2 



£ ( €-ij ki Tj / 



one can rewrite the statistical mixing Eq. (61) as 

2 



i 



7T 2 ft 2 e CT iY(cr) 



s=l Ln s =l J 

To obtain an explicit expression for N(cr), one uses the following property of the Laguerre 
polynomials, 



E^)A' = ?r ^p I 



exp 



A 



A - 1 



through which, for e CTs < 1, with A = — e CTs , one obtains 

oo 1 

EH"-)"' t[<" + 2(-i)'£] = «+^: 



exp 



n s =l 

which results into 



|A| < 1; 



e 2 + 2(-l) s £ 
1 + 



(65) 



(66) 



1 



exp 



£ 2 sinh (cr) + 2£ sinh (ecr) 
cosh (a) + cosh (ecr) 



(67) 



2 7r 2 ft 2 iV(cr) [cosh (a) + cosh (ecr)] 
The normalization integral is evaluated over the infinite four-parameter phase-space, 
{{ r i> { r 2, and hence one obtains 



N(a) 



1 



2 cosh (cr) — cosh (ecr) 



such that 



1 cosh (cr) — cosh (ecr) 
7r 2 ft 2 cosh (a) + cosh (ecr) 



exp 



£ 2 sinh (a) + 2£ sinh (ecr) 
cosh (cr) + cosh (ecr) 



(68) 



• (69) 
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Since the above result depends on time-invariant quantities, £ 2 and C (c. f. Eqs. (39)), 
one can obtain the corresponding probability distribution relative to coordinates Qi(<-> ri) 
and Q 2 (^- r 2) by integrating the above thermalized Wigner function over k\ and k 2 , as 

P{<r,r u r 2 ) = dkj dk 2 pJKa; £ 2 , £). (70) 

'' — oo J — oo 

The distortion over the probability distribution resulting from NC effects can be obtained 
by subtracting "p( e_s> °) from V. The result follows from the axial invariance (V(cr; r\, r 2 ) = 
V(cr; Q 2 )) which can be visualized in the Qi — Q 2 plane, as depicted from Fig. [6] The results 
are described in terms of an increasing scale of the thermodynamic parameter o = Ml/ (ksT). 
For illustration purposes, the pictures are exemplified for e = j/fl equals to 0.1 and 0.5. 

To provide a more accurate estimate of the contribution of noncommutativity on the 
missing information, one examines the quantum entropy and the quantum mutual infor- 



mation of the system described above (c. f. the expression for S\ 2 , Eq. (58)). The linear 
entropy, Si 2 (a), and the mutual information, Ii 2 (a), as well as the corresponding missing 
information described in terms of AS12 ~ 5*12 — S^ ^ and A/ 12 ~ 7 12 — , are depicted 
in Fig. [7j 

One notices that the isotropic 2D harmonic oscillator does exhibit an effect of noncom- 



mutativity through the statistical mixing, Eq. (61). 

At high temperatures, the corrections due to noncommutativity reflect themselves into 
the missing information given by A £12 oc e 2 cr 2 and by A/12 oc e 2 a, which are therefore small. 
As the parameter a decreases, the system can access quantum states with larger quantum 
numbers, the linear entropy approaches the unity, and the system exhibits the behavior 
of a maximal statistical mixing, i. e. the classical limit. For extremely low temperatures, 
ksT Ml (a 1), the system accesses predominantly quantum states with small quantum 
numbers. This corresponds to a decreasing level of mixing and hence to quantum behavior, 
despite NC contributions. For e < 1 (7 < O) the usual pattern of missing information is 
maximally modified by the NC element at intermediate scales (1 < & max < 4), as one can 
notice from numerical results depicted in Fig. [8} There is a transition regime where the NC 
effects are maximized. 

One can also see from Fig. [6] that the distortion of the probability distribution and an 
increase in the of missing information, AS12 and A/12, coincide. The classical limit (cr 1) 
and the harmonic oscillator ground state (cr ^ 1) are obtained at the asymptotic limits of T 
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and a, respectively. Any additional decoherence effect that is present in the evolution of a 
quantum superposition of state vectors with well-defined quantum numbers are therefore a 
natural extension/consequence of the orbital distortion from Fig. [6j which reflects an overall 
increase of the entropy of the system. 



A. The thermodynamic limit 

A large system of 2D NC harmonic oscillators which is in thermal contact with an environ- 
ment, with temperature T, corresponds to a canonical ensemble. The microstates occupied 
by the system are described by n , where E ni>n2 denotes the (star-gen) energy of the system 
in a given microstate. These microstates can be regarded as a system of discrete quantum 
states. 

The classical partition function is obtained when the trace is expressed in terms of the 
coherent state vector, p^. n , and when quantum uncertainties are negligible. Formally, one 



has 



Z(a) 



|-oo „ + oo 



+ 00 „ + oo 

&T\ \ dr 2 I dki I dk 2 exp 



00 



exp 



"1,712 =0 

e - ' 



E, 



ni,n 2 



k B T 



TTW 

k B T 

00 



P W 



exp 



[— cr] 2J exp [— cr(l + e)ni — cr(l — e)n 2 ] 



n\,n2=0 



1 



1 _ e -<T(l+e) I _ e -(r(l-e) 

1 1 

2 cosh (cr) — cosh (ecr) 



(7 



(71) 



which is consistent with Eq. (68). 



The partition function, Z(o~), allows for obtaining the thermodynamic variables of the 
system, such as the internal energy: 



U(a) = -hn-^-\n [Z(a)]; 



(72) 



the Boltzmannn entropy: 



S k (a) = -k B {In [Z{a)\ - aU(a)/hn} ; 



(73) 



and the heat capacity: 



r%2 

C v = k B a 2 — \n[Z{a)] 
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(74) 



By substituting the result from Eq. (71) into the above definitions, one obtains 

U(cr) — M] Sinh ^ _ 6 Sinl1 ^ e<J ^ (75) 
cosh (a) — cosh (ecr) ' 

(cosh (cr) — cosh (ecr)) In [2 cosh (a) — 2 cosh (ecr)] + a [sinh (cr) — sinh (ecr)] 

^fe(cr) = —KB 



cosh (cr) — cosh (ecr 
and 



(76) 



^ _ ^ cr 2 (1 + e 2 ) [cosh (cr) cosh (ecr) - 1] - 2cr 2 e sinh (cr) sinh (ecr) 

(cosh (cr) — cosh (ecr)) 2 

from which the NC corrections are evinced through the explicit dependence on e = 7/O. 

The NC effects on the above defined thermodynamic variables, Sk and C v , are depicted 
in Fig. [9j The qualitative behavior of the Boltzmannn entropy, Sk, is consistent with that 
of the quantum linear entropy, in the sense that the NC features add up to the entropy 
The heat capacity, C v , shows an anomalous dependence on the thermodynamic parameter 
cr such that the NC corrections may increase or decrease with respect to C^ - ^ (c. f. the 



second plot of Fig. |9J). For these thermodynamic variables, the effect of noncommutativity 
is maximal for intermediate scales of cr. This offers a novel possibility for measuring the NC 
effects at low temperatures. 



V. AXIAL SYMMETRY AND THE r?-ZEEMAN EFFECT 



It can be speculated whether the NC effects exhibited by the isotropic harmonic oscillator 
can be factorized from the Wigner distribution if one chooses the appropriate representa- 
tion for the quantum operators. Presumably, the axial symmetry of the problem might be 
relevant to obtain a suitable representation. 



From the point of view of the quantum system described by the Hamitonian, Eq. (32), the 



noncommutativity can be re-interpreted as an analogous of the Zeeman effect for a charged 
rotating particle in the presence of an external magnetic field B. 

In this context, it is helpful to consider the three-dimensional (3D) extension of the 
previously discussed NC problem. 

The Hamiltonian for the problem can be written as 



#ffo(q,p) 



P , 1 22 

— moo q , 

2m 2 H 



1 

2m 



+ 



muj 



Q 



2h 



(o x n) 



(78) 
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since we have used a 3D generalization of the SW map given by 

Qi = Qi - ^i jk 9jtl k , pi = fli + ^ e ijkVj Qk , (79) 
that reflects the following NC rules, 

[qi, qj] = ie ijk e k , [q^Pj] = iS^h ( 1 - ^-^ J , = le^k , z,j = 1,2,3 (80) 



4n 2 f 

where we have set p = A = 1, is the Levi-Civita tensor of rank three, 6 = {0i, 62,03), 
and r] = (771, T]2,T]3). The 2D version of the problem can be recovered from above relations 
by setting 9 1%2 = Vi,2 = 0, 3 = 0, and r] 3 = rj. 

Even for a free particle, NC effects on the Weyl-Wigner Hamiltonian, H^ Q , will be present 

in the particle dynamics through its kinetic energy: 

.2 



1 12 



One can use that 

77 ■ ( Q x II) = 77 • L = t]L cos (if) = r]L z , (82) 

and 

(r, x Qf = rfR 2 , (83) 

with R = I Q\ sin (</?) = \/ Qf + Q$, where ip = arctan Rj Q3 is the angle measured on the 
plane of axial symmetry along z, defined by an arbitrary direction of r/ (that in this case 



corresponds to 77 = rjz). The axial symmetry ensures that the Hamiltonian Eq. (81) can be 
re-written as 



Dropping the q-dependent potential, the Moyal siar-product properties lead to the 
stargenvalue equation written as 



H H o(p)*^q)=H v (p)* tl * h ^q) = H^(U,Q)-kn^(Q)=E^(Q). (85) 

The second-order correction due to the ?7 2 -term has to be considered in order to set a complete 
map of the Zeeman effect [52]. The replacement of the parameter rj by the coefficient 
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eBh/c into Eq. (84) maps into the dynamics of a particle with charge, e, moving in a plane 



perpendicular to the magnetic field B = Bz. 

The eigenfunctions are obtained in terms of cylindrical coordinates (R,(p,z) as ^(Q) = 
ipK,i(R, f) Z{z). Since Z{z) cx exp (—ik z z/h) is not relevant for the present purpose, one has 



i^(II,QK*(Q) = lH?(R,<p)- 



2m dz 2 



if> K AR,(p) z 



E K/ + ^ ) ^(R, <p) Z. 



The remaining eigenvalue equation becomes 

h 2 ( d 2 



2m \dR 2 ' RdR R 2 dip 2 J 2m dip 8mh 2 
E K) £ ip Kt i(R,ip), 



(86) 

^K,t{R,w) 

(87) 



where we have identified IL, = —ih(d/dQi), with i — 1,2, and L z = —ih(d/dip), and the 
simultaneous eigenfunction of and L z , i[) K AR,<f), can be obtained by separation of 
variables as ip K ,t(R, <f) = f*A R ) e%h '• 

The resulting ordinary differential equation for f K ^(R) is given by 

h 2 ( d 2 id ^ °- - 2 



\ It] vr E 

2m V OR 2 'RdR R 2 2m 8mh 2 K ' 1 



for which, after some mathematical manipulations, we obtain the solution 

Mi 

2 

exp 



V „ 2 



and the energy spectrum depending on the quantum numbers, k and £: 



E K 



V 

2m 



(2K + \e\-e + i). 



(88) 



(89) 



(90) 



If one considers the harmonic oscillator potential at the NC Hamiltonian from Eq. ( 78 ) 



the result of the above analysis would be valid up to second-order in 9 and 77, and would 
lead to the energy spectrum: 



E K>e = 2ha(3(2K + \£\ + 1) - h-f£, 

which relates to the solution described in the previous sections through 2k + \£\ 
and £ — n-i — n\. The normalized eigenfunction, in that case, would be given by 



(91) 



ijj K AR,<f) 



a 

hp 



R 2 L W I " # 



a 



exp {it ip). 



(92) 
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with a and f3 redefined up to second-order in 9 and r\. The relevance of separating the 
solution in cylindrical coordinates is that the re-defined quantum numbers introduce energy 
eigenfunctions and eigenvalues with well-stablished physical meaning. To sum up, the ana- 
lytic solution of the problem that contains second-order corrections in B (or, equivalently, 
7] 2 , in the NC version) is provided in terms of associated Laguerre functions L«'. Notice that 
the solution f K (R) has been presented in terms of a regular power series expansion (see, for 
instance, the resolution of problem 11.6, pag. 632, at Ref. [52]). Thus, this result is useful 
even for the ordinary QM problem. 

VI. CONCLUSIONS 

In this work we have examined the implications of noncommutativity on quantum beat- 
ing, scale dependent decoherence and dissipation effects related to the loss of information 
of quantum systems. Considering the time-evolution of the 2D harmonic oscillator in the 
NC phase-space through the Groenwold-Moyal description of the Wigner function, we have 
reconstructed a suitable framework for investigating NC effects. It is shown that the har- 
monic oscillator on the NC plane approaches the classical limit, and exhibits well marked 
quantum effects like state-swapping, quantum beating, and loss of quantum coherence to 
some extend. 

We have shown that phase-space noncommutativity may induce the destruction of the 
original features of the wave function for the commutative problem. The NC terms distort 
the 2D harmonic oscillator quantum orbitals through modifications onto their corresponding 
probability distributions. To study these features, the missing information quantified by the 
linear entropy was computed in the NC context. The linear entropy and the quantum 
mutual information of the 2D NC harmonic oscillator were examined and the conditions for 
identifying the maximal missing information due to the phase-space noncommutativity was 
obtained. 

Finally, we have also investigated the thermodynamic limit of the quantum system by 
discussing the behavior of the Boltzmannn entropy and of the heat capacity derived from 
the partition function obtained for the thermalized 2D NC harmonic oscillator. A corre- 
spondence between the NC effects over quantum and classical variables were evinced, and 
it was shown that the most suitable range of temperatures for detecting noncommutativity 
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is around ~ hVt/kB (c. f. Eq. (37)). 

As a complementary issue, the 2D NC harmonic oscillator problem was rewritten in 
order to depict an axial symmetry so to ensure a suitable factorization of the NC effects. An 
analogy with Zeeman-effect has been drawn, resulting in a description of the system with 
well-defined NC quantum number and a factorable wave-function. The results allow us to 
decouple the non commutative effects such that it resembles a Zeeman-like effect. 

One additional feature of the axial symmetry of the problem deserves some comments. 
The introduction of an intrinsic anisotropy into the quadratic potential of the NC harmonic 
oscillator would not qualitatively change the obtained results. Besides suppressing the pos- 
sibility of decoupling the NC effects that break the isotropy through the introduction of 
different elastic constants, k± oc ojf and ki oc w|, would lead to different quantum beat- 
ing patterns. In this case, in order to recover the initial configuration of the phase-space, 
the system would require a number of cycles proportional to the minimal common integer 
multiple of oscillation frequencies fix, Q 2 , 7i and 72 . 

One could also consider the possibility of introducing a linear perturbation to the har- 
monic oscillator in the coordinate variables (Qi or Q2) in order to produce some kind of 
anisotropy between Qi and Actually, this would only produce a relative shift in the 
phase-spaces ((Qi,H2) or (Q2,Hi)) likewise the effect of gravity in the unidimensional har- 
monic oscillator in z direction. 

To summarize, our results suggest that noncomutativity effects can be interestingly con- 
sidered when addressing the issues of interface between quantum and classical descriptions 
of nature [5?]. As it happens, for quantum superpositions, that are irreversibly affected by 
their surroundings, the quantum collapse due to noncommutativity may reflects the loss of 
coherence of pure states that are continuously transformed into a statistical mixture. In the 
general context of the discussion of mechanisms to induce decoherence and loss of informa- 
tion [58J , noncommutativity is shown to play a relevant role in inducing the transition from 
quantum to classical behavior. 
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FIG. 1: Time evolution of the phase-space coordinates, (Qi(i), III (t)) and (Q2(t),H2(t)), for the NC har- 
monic oscillator. The first plot line refers to the phase-space elliptical orbits of commutative harmonic 
oscillators (as if one had set e = in the NC map). From the second to the forth plot lines one has set 
arbitrary values for e, e = 1/4, 1/10, and 1/100 respectively. The precession motion describes close orbits 
for these cases. One has used a BlueGreenYellow scale in order to denote r = (blue) and r = 2?r/fi 
(yellow), such that orbits start and finish at (x, tt x ) = (y, ir y ) = (1/V2, 0). The last pair of plots corresponds 
to e = 1/27T, which does not result into a closed orbit since it is not a rational number. By convenience, one 
has set a = f3, that is equivalent to muj = h = 1 
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FIG. 2: Quantum beating for the NC harmonic oscillator described by the traced-out state vector 
Pn},n y (Qi, H; f) in the Qi— IIi plane, with n y = 2 and n x = 1 (first column), 2 (second column) and 
5 (third column). At time r = the Wigner function is assumed to be centered at the origin (0, 0). 
One has considered time intervals such that r = fc7r(8eu;) -1 , with k = 0, 1, 2, 3 and 4 in order to 
reproduce its time evolution through the sequence of five plots. Notice that the prepared quantum 
state in Q\ — Yi\ is projected into the traced-out quantum state, by reproducing its wave pattern at 
r = 7r/(2ew). It corresponds to the quantum beating effect with frequency Ubeat = 2j = 2eoo. The 
contour plot follows a rainbow scale (from red which corresponds to 1, to violet which corresponds 
to 0). 
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Quantum beating for the NC harmonic oscillator described by the state vector 
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FIG. 3: 

column), 2 (second column) and 5 (third column). The resulting Wigner function in the Qi — Hi 
plane is supposed to be centered at (<3i,IIi) = (1, 0) at time r = 0. One has considered time 
intervals such that r = /c-7r(8eu;) -1 , with k = 0, 1, 2, 3 and 4. Once again one notices a quantum 
beating effect with frequency ujbeat = 27 = 2eu. Beside the analogous effects described in Fig. [2j 
the NC parameter e introduces the precession motion that follows the phase-space maps like those 
ones in Fig. [T] 
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FIG. 4: Mutual information, I\2, for the NC harmonic oscillator described by the state vector 
for several values of the quantum numbers n x and n y . 
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5: Approximated decoherence profile for tiny values of the NC parameter e, depicted for 
increasing values of the mutual information, I\2- 
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FIG. 6: Probability distribution (quantum orbital) distortion for the NC harmonic oscillator state 
vectors with radial symmetry in the Q\ — Qi plane. For higher temperatures, i. e. for fcgT S> hco 
(a <C 1), the NC effects are erased, and the quantum behavior of the thermalized state vector, p^, 
prevails over the tiny NC quantum corrections. For very low temperatures, i. e. for ksT <C foj 
(a S> 1), the wave function collapses into the classical limit. The usual pattern of the quantum 
orbitals is maximally modified by the NC element at intermediate scales, namely 1 < (Tmax i$ 4. 
The plots correspond to e = 0.1 (l s *-column) and 0.5 (2 nd -column) with a assuming the values of 
0.1, 0.2, 0.5, 1, 2, 5 and 10 (decreasing temperature scale). 
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FIG. 7: Linear entropy (black lines), S^c), and mutual information (blue lines), 112(0"), and 
corresponding missing information given in terms of A5i2(cr) and AIi2(cr), for the isotropic 2D 
NC harmonic oscillator. For higher temperatures, ksT ^> hj (a <C 1), the NC effects are not 
effective. For very lower temperatures, fc^T C % (ff > 1), the wave function collapses into the 
classical limit. The maximal values of A1S12 (cr) and A/i2(<r) are obtained at intermediate scales of 
a (see the correspondence in Fig. [8]). One has considered e = (solid lines), 0.1 (dotted lines) and 
0.5 (dashed lines). 
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FIG. 8: Numerical solution for the a parameter that maximizes the missing information (from 
Fig. [7]) as function of the NC parameter e = 7/fi. The results are for the maximal values of 
ASi2(c) (black lines), and Al^cr) (blue lines). 
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FIG. 9: Thermodynamic variables: Boltzmann entropy (black lines), Sk(cr)/kB, and heat capacity 
(blue lines), C v (a)/kB, for the statistical mixing pf h . It reproduces the thermodynamic limits of the 
isotropic 2D NC harmonic oscillator quantum behavior. The distortion due to the NC effects are 
quantified through ASk(cr)/kB and AC v (a)/kB from the second plot. One notices that C v (a)/kB 
exhibits an anomalous behavior that inverts the sign of AC v (a)/kB at intermediate scales of a. 
The plots are for e = (solid lines), 0.1 (dotted lines) and 0.5 (dashed lines). 
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